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Abstract 


A set of mathematical formulation is adopted to study vapor deposition from 
source materials driven by heat transfer process under normal and oblique 
directions of gravitational acceleration with extremely low pressure environment 
of 10“ 2 mm Hg. A series of time animation of the initiation and development of 
flow and temperature profiles during the course of vapor deposition has been 
obtained through the numerical computation. Computations show that the process 
of vapor deposition has been accomplished by the transfer of vapor through a 
fairly complicated flow pattern of recirculation under normal direction 
gravitational acceleration. It is obvious that there is no way to produce a 
homogeneous thin crystalline films with fine grains under such a complicated 
flow pattern of recirculation with a non-uniform temperature distribution under 
normal direction gravitational acceleration. There is no vapor deposition due 
to a stably stratified medium without convection for reverse normal direction 
gravitational acceleration. Vapor deposition under oblique direction 
gravitational acceleration introduces a reduced gravitational acceleration in 
vertical direction which is favorable to produce a homogenous thin crystalline 
films. However, oblique direction gravitational acceleration also induces an 
unfavorable gravitational acceleration along horizontal direction which is 
responsible to initiate a complicated flow pattern of recirculation. In other 
words, it is necessary to carry out vapor deposition under a reduced gravity in 
the future space shuttle experiments with extremely low pressure environment to 
process vapor deposition with a homogenous crystalline films with fine graims . 
Fluid mechanics simulation can be used as a tool to suggest most optimistic way 
of experiment with best setup to achieve the goal of processing best nonlinear 
optical materials. 
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I . Introduction 


In recent years, a great deal of interest has been directed toward the use 
of organic materials in the development of higher efficiency optoelectronic 
devices. There is a myriad of possibilities among organics which allows a great 
deal of flexibility in the design of unique structures with a variety of 
functional objectives. The use of nonlinear optical (NLO) organic materials as 
thin film waveguides allows full exploitation of their desirable qualities by 
permitting long interaction lengths and large power densities with modest power 

input(Nayar and Winter, 1990). There are several methods in use to prepare thin 

♦ 

films, such as Langmuir -Blodgett (Carter et al. , 1988; Kajzar et al. , 1983; Kajzar 
and Meissier, 1988), growth from sheared solution or melt(Thakur and Meyler , 
1985; Thakur et al . , 1986; Seymour et al., 1985), and melt growth between glass 
plates (Ledoux et al . , 1986). This method, for preparing thin films, has shown 
some evidence that microgravity processing with extremely low pressure of 10 2 
mm Hg may be beneficial. 

Epitaxial growth on ordered organic and inorganic substrates and variations 
in processing conditions have been used to prepare highly oriented 
polydiacetylene and phthalocyanine films. The relationship of processing 
conditions to uniformity in thickness, degree of orientation, and nonlinear 
optical properties are very important. Also, improvements in fabrication 
techniques for growth of thin films by vapor deposition would be useful. Time 
variation of the initiation and formation of flow field and temperature 
distribution during the vapor deposition process under extremely low pressure 
environment can provide a better understanding in the framework to determine 
whether nonlinear optical thin films should be grown in space to improve 
uniformity in thickness and homogeneity if ground-based research and theoretical 
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studies indicate that this is feasible. 


NASA Marshall Space Flight Center is currently engaged in the research of 
using a vapor deposition method to process nonlinear optical materials by growing 
thin crystalline films of diacetylene monomer (Paley et al., 1992). A better 
understanding of flow profiles and temperature distribution for vapor, which is 
produced from source material and propagates through the major driving mechanisms 
of heat transfer with natural convection and the diffusion process, can greatly 
enhance and improve the technique of nonlinear optical material processing. 

In this study, an initiation and development of time - dependent , three- 
dimensional fluctuations of flow profiles and temperature distribution driven by 
heat transfer with a natural convection under normal and oblique directions of 
gravitational acceleration with extremely low pressure environment have been 
considered. Purpose of this study is to investigate the probability of 
processing vapor deposition under normal gravity with various directions to 
produce a homogenous thin crystalline films with fine grains . Results of present 
study will provide the critical guidelines for future space shuttle experiments. 
Size of container adopts in this simulation will be served directly to facilitate 
the infrastructures to be used for both ground-based and space shuttle 
experiments . 

II. Mathematical Formulation 

Governing equations, initial and boundary conditions for the present study 
of convection- driven material processing are illustrated. 

(II -A) Governing Equations 

Consider a co-axial cylinder with an inner radius of and an outer radius 
of R 2 . The heights of the outer and inner cylinders are L 2 , and (L 2 - Li) , 
respectively. Figure 1 shows the geometrical configuration of the co-axial 
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cylinder. For the purpose of considering the mathematical formulation, 
cylindrical coordinates (r, 9, z) with a corresponding velocity components (u, 
v, w) , and corresponding components of gravity acceleration (g r , g$, g z ) i£ 
adopted. Figures 1(A) and 1(B) show the geometrical configuration of co-axial 
cylinder in the r-z plane and r -9 plane, respectively. There are two temperature 
zones in this co-axial cylinder. The external walls of the outer cylinder are 
surrounded by water jacket which maintains a constant temperature of 120° C while 
the walls of the inner cylinder keeps a constant temperature of 30° C. A low 
pressure air with the pressure of 10~ 2 mm H g (= 1.333 N/m 2 ) is filled in the gap 
between the outer and inner walls of the co-axial cylinder. The flow profiles 
and temperature distributions of the low pressure air are mainly governed by the 
natural convection of the flow pattern in this problem. The governing equations 
can be shown as follows : 

(a) Continuity Equation 


l_d_ 
r dr 


( ru ) 


+ 


1 dv 

r ae 



(i) 


(b) Momentum Equations 


. du du 
p< at +u 3? 




2 du 

+w 


dp 


+ ,(VWL-4-g> 


( 2 ) 


_ / dv 

p( Tc 


+ U 


dv 

dr 


v dv 

r ae 


HZ + ^.fZ) 

r dz 


1 dp 

r ae 




3 



(4) 




=-||- p sr 2 [ 1 -P ( r-Tj]+nV 2 w 


where 




dT 


1 32 + . & 


r dr r 2 d0 2 3z 2 


and 


4 


(g Z '&Q'9z'> = flr(sini|rcos0, -sim|rsin0 # -cosi|r) 


(c) Energy Equation 


*“f ‘Tf *'I> ' 


(5) 


Here, p, p, T, /i, k, and g denote density, pressure, temperature, viscous 
coefficient, heat conduction coefficient of air, and background gravity 
acceleration, respectively. tf> is the angle between the gravity vector and the 
axis of cylindrical cylinder. In this case, To, and p m stand both air temperature 
and density at the conditions of low pressure at 10" 2 mm Hg and heat sink at the 
wall of the inner cylinder at T = 30° C. In these mathematical expressions, 
Boussinesq approximation is employed to formulate the relationship between local 
fluctuations of density and temperature through the thermal expansion coefficient 
which replaces the equation of state. 

(II-B) Initial and Boundary Conditions 

The initial and boundary conditions have been assigned explicitly to solve 
governing equations (11) to (5) simultaneously. 

(a) Initial Condition 

At time t = 0, all the flow profiles vanish, i.e. 

u = v = w = 0 everywhere (&> 

The temperature distribution of fluid inside the container is 30°C, and 
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outside the container is expressed by 


T = 30° C at r < R x , 0 < 9 < 2tt, (L 2 - L x ) < z < L 2 (7) 

T = 120° C at z < 0; r > 0; 0 < 9 < 2n (8) 

atz>L 2 ;r>O;O<0<27r (9) 

at 0 < z < L 2 ; r > R 2 ; 0 < 9 < 2tt (10) 

Pressure profile inside the container is 

p — 10" 2 mm Hg = 1.333 N/m 2 (11) 

(b) Boundary Conditions 


At time t > 0, temperature distribution shown in Equations (7) to (10) , and 
pressure profile shown in Equation (11) remain to be true. 

For flow profiles, no penetration and no slip conditions along the surfaces 


of solid wall shall apply, namely 

u = v = w = 0 (12) 

along the surfaces of r = R x , 0 < 9 < 2n , (L 2 - L x ) < z < L 2 ; (13) 

r > 0, 0 < $ < 2tt, z — 0 ; (14) 

r>0, 0 < 9 < 2ir , z = L 2 ; and (15) 

r - R 2 , 0 < 0 < 27 t, 0 < z < L 2 (16) 


A co-axial cylindrical container designed to be used directly in the future 
space shuttle experiments, illustrated in Figure 1, with following specifications 
are assigned: L x =6 cm, L 2 =12 cm, R x =0.9 cm, and R 2 =l . 5 cm. Following physical 

parameters were used: gravity acceleration g Q = 9.81 m/s 2 , air pressure=l . 333 
N/m 2 , air density=l . 71xl0" 5 kg/m 3 , air viscosity coef f icient=2 . 0x10 5 m 2 /s, 

coefficient of thermal expansion for air (/3)=2 . 87xl0" 3 1/K and Prandtl number for 
air=0 . 72 . 

to investigate the most significant driving parameters in this study, let 
us consider the inertia, buoyancy and viscous terms, 

w»w/L - g z ^AT - vW/8j 2 (17) 
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respectively, which are on the same order of magnitude in momentum equation shown 
in Equation (4); and convection and conduction terms, 

w/LAT ~ a AT/$ t 2 (18) , 

respectively, which are also on the same order of magnitude in energy equation 
shown in Equation (5). Here, S x is the thermal boundary layer thickness; L is 
the height of container along z-axis; u = p/p is the kinematic viscosity; and a 
= k/pCp is the thermal diffusivity. From convection and conduction terms in 
Equation (18) ; from inertia and buoyancy terms in Equation (17) ; and from 
bouyancy and viscous terms in Equation (17) , we can obtain induced convective 
velocity of flows as follows: 

w ~ qL/ 6 t 2 (19) 

w - (g z /3LAT) 1/2 (20) 

and w ~ (g z ^6 T 2 AT)/i/ (21) 

respectively. By equating Equations (19) and (21) , we have the relationship of 
S t /L ~ Ra" 1 '* (22) 

based on the balance between bouyancy and viscous forces. By equating Equations 
(19) and 20) , we have another relationship of 

8 T / L ~ (Ra- Pr)” 1/A ~ (Gr)" 1/A (Pr) _1/2 (23) 

based on the balance between buoyancy and inertia forces. Here Reyleigh number 
(Ra) , Grashof number (Gr) , and Prandtl number (Pr) are defined as 



Ra - g z ^L 3 AT / av 

(24) 


Gr - g z /9L 3 AT/y 2 - Ra/Pr 

(25) 

and 




Pr = v/cl 

( 26 ) 

Since 

the heat transfer coefficient scales as k/6 T (Bejan, 

1995) , the Nussel 

number 

(Nu) varies as 



Nu = hL/k - Ra 1M 

(27) 
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from Equation (22) based on the balance between buoyancy and viscous forces; and 
Nu varies as 

Nu - (Ra- Pr) 1/4 ~ (Gr) 1/4 (Pr) 1/2 (28) ^ 
As indicated by Bejan (1995), both theoretical and laboratory measurement show 
that Equation (27) is particularly applicable for Pr > 1 fluids, while Equation 
(28) is dominated by Pr « 1 fluids. For the case of present research, the 
dominated parameter is Rayleigh number rather than Grashof number (Bejan, 1995) . 

In this study, it is shown later in the numerical simulation that the 
maximum induced flow velocity driven by convection is 0.2 m/s. Reynolds number 
(Re = wD/i/ where D is the diameter of container) for the flow pattern based on 
this maximum flow velocity is 150 which is far below the possibility to produce 
turbulent flow fields. Grashof number (Gr) and Rayleigh number (Ra) for the 
present case study are 8 X 10 6 and 7.8 X 10 6 , respectively, under normal gravity 
condition. It is well-known that the condition for the onset of convection with 
fluid layers heated from below is Ra > 1.1 X 10 3 (Bejan, 1995) . In other words, 
the total elimination of the onset of induced convection with the proposed size 
of container is possible only to perform the space experiment of material 
processing under the reduced gravity environment of 10' 4 g 0 and below. 

III. Mathematical Simulation of Flow Profiles and Temperature Distributions 

Equations (1) to (5) subject to initial conditions shown in Equations (8) 
to (11) , and boundary conditions shown in Equations (12) to (16) , have been used 
to solve time -dependent, three-dimensional flow profiles and temperature 
distributions numerically. The computer algorithms employed in this study have 
been developed and illustrated in our earlier studies (Hung et al . , 1992 a, b, 
c, d; 1993 a, b; 1994 a, b, c, d, e) , and will not be repeated in this paper. 
The time step is determined automatically based on the size of grid points and 
the velocity of flow fields. As the thickness of boundary layer is inversely 
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proportional to the square root of Reynolds number, size of grid points shall be 
adjusted according to the flow fields of Reynolds number. In other words, the 
size of grid points shall be smaller as the flow fields of Reynolds number 
becomes greater. Fail of adjusting the size of grid points based on the 
increasing of Reynolds number will result in the numerical instability to the 
extreme severity of the local oscillations ( Kitchen et al., 1978; Kitchens, 
1980) . Special care shall be attempted to eliminate these numerical 
instabilities by properly adjusting the size of grid points which eventually 
contribute to the adjustment of time steps. It is needless to say that the 
development of a real physical instability will certainly occur at the higher 
Reynolds number because of less viscosity effect contributing to the damping of 
flow disturbances. It is utterly important for us to distinguish the differences 
between numerical and real physical instabilities, particularly for the flows 
with high Reynolds number. 

In this study, two sets of the size of grids, with the grids of 17 X 34 X 
42 and 23 X 34 X 82 along (r,0,z) coordinates, have been designed to carry out 
numerical computations . As the Reynolds number of 150 with maximum induced 
velocity of 0.2 m/s for the present case, there is no difference concluded for 
numerical results obtained for either grids chosen. By using the grids of 17 X 
34 X 42 in the computation of flows with higher Reynolds number, such as Re > 
1500 with induced velocity of 2.0 m/s and higher. An extremely severe numerical 
instabilities of flow field is resulted. 

With grids of 17x34x42, the average time step determined from the size of 
grid points and velocity of flow fields is 0.0007 s, while the CPU time required 
by the Supercomputer CRAY II to execute one time step is 3.019 s. The total CPU 
time consumed for each of the three cases in this study is 12000 s CRAY XI 
computer time. 
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In this study, direction of background gravity along three different 
directions with tf> = 0° , 45° and 180° along the r-z plane with 0=0° are chosen 
as the examples to investigate the characteristics of convection- driven flow and 
temperature profiles. Results of present study will be used as the guidelines 
to facilitate the space experiments to be carried out with the container 
specified in this paper. 

Various references have been cited by Bejan (1995) in the study of induced 
convective heat transfer with fluid layers heated from below, top, side, oblique, 
etc., under the normal atmospheric pressure. As mentioned in the Introduction, 
purpose of this study is to investigate time variation of the initiation and 
formation of flow field and temperature distribution during the vapor deposition 
process to determine the possibility of processing nonlinear optical thin films 
under the extremely low pressure environment of 10~ 2 mm Hg. Study of induced 
convective heat transfer under such a low pressure environment has not been done 
earlier, and there is no available literature for comparison. However, there is 
a qualitative agreement of both flow profiles and temperature distribution for 
the induced convections under normal atmosphere pressure and 10 2 mm Hg low 
pressure . 

(A) Convection Driven by Normal Direction 0°) Gravitational Acceleration 
Time animation of the development of flow fields with normal direction (^>= 
0°) gravitational acceleration is illustrated in Figure 2 for convection- induced 
flow vectors (in terms of directions and magnitudes of convection induced- flow 
velocities) in r-z plane. The maximum velocity for each subfigure is specified 
separately at the corresponding subfigure. For the convenience of comparison, 8 
subfigures with time sequence at t= 0.22, 0.58, 1.29, 1.72, 1.99, 2.39, 2.70, and 
3.10 s throughout this paper are illustrated. Following results are concluded 
in these series of subfigures in Figure 2: (a) A convection induced downward flow 
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toward the bottom of outer cylinder with inflow stream is initiated at the bottom 
of central column at time t<0.2 s (not shown), (b) A convection induced downward 
flow along the surface of cylindrical central column with inflow from the surface 
of outer cylinder toward the surface of cylindrical central column (at upper 
section of container within the slender gap between co-axial cylinders) is also 
initiated during the time period of 0.2 s<t<0.4 s. (c) A convection induced flow 
patterns of upward flow along the surface of outer cylinder and downward flow 
along the surface of inner cylinder (central column) are gradually formed (at 
upper section of container within the slender gap between co-axial cylinders) 
with intensified flow velocities accompanied by a stream of inflow from the 
surface of outer cylinder to inner cylinder during the time period of 0.2 s<t<0.9 
s. (d) A tendency of convection induced inflow pattern toward center of axial 
coordinate at upper section of container with heights in z>2.7 cm, and outflow 
pattern toward the surface of outer cylinder at lower section of container with 
heights in z<2 . 7 cm are gradually formed for time t>0.9 s. (e) It shows a 
convection induced flow pattern of two recirculation flows with centers at 
(r c , z c ) = (±l . 9 , 2.7) cm where r c and z c denote the centers of recirculation flow 
along r- and z- coordinates, respectively. Direction of left-hand-side 
convection induced recirculation flow is clockwise while that of the right-hand- 
side convection induced recirculation flow is counter-clockwise circulation for 
time t>l .3 s . 

For the purpose to illustrate the flow patterns of convection induced 
inflow and outflow, Figures 3 and 4 show time animation of initiation and 
development of convection induced inflow in the region below the central column 
(at z=5.55 cm), and convection induced outflow in the region near the bottom of 
outer cylinder (at z=0.45 cm), respectively. Figures 3 and 4 also show 
convection induced flow vectors (in terms of directions and magnitudes of induced 
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flow velocities) in r -9 (x-y) plane. The maximum velocity for each subfigure is 
also specified separately at the corresponding subfigure. It is interesting to 
emphasize that there is an induced inflow region with flow from the surface of 
outer cylinder inward toward the central direction with location at the upper 
section of container with heights in z>2 . 7 cm. It is also an induced outflow 
region with flow from the center of axial coordinate outward toward the surface 
of outer cylinder at the lower section of container with heights in z<2.7 cm. 
Before the induced flow is developed to such a clear-cut regions of upper part 
with inflow and lowerpart with outflow for the flow approaching to steady- state 
case, there is an induced transient flow pattern at the beginning of flow 
formation. Figure 3 shows that there is an induced temporary outflow formed at 
the exterior region of container near the surface of outer cylinder at the very 
beginning of flow initiation. This induced temporary flow region disappears 
quickly as the flow development near the steady-state. Obviously, the major 
vapor deposition at the bottom surface of the central column are contributed by 
the induced inflow from exterior flow zones which were brought upward along the 
surface of outer cylinder from the evaporation of source materials stored at the 
bottom surface of outer cylinder. In other words, the flow analysis indicates 
that the vapor deposition at the bottom surface of central column was not a 
result of direct straight upward diffusion from the evaporation of source 
materials stored at the bottom surface of outer cylinder. It is rather a vapor 
deposition process resulted from a fairly complicated induced flow pattern of 
recirculation. 

Figure 5 shows a time animation of the initiation and formation of 
temperature profiles due to induced convection in r-z plane at 0=0° and 180° 
during the course of the entire induced flow region development. The curves show 
in Figure 5 are the isothermal lines with values indicated with temperature in 
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°C at the corresponding locations. Comparisons of Figures 2 to 5 can be 
concluded as follows: (a) At time t-o (not shown), there is a very sharp 
temperature varying from 30 to 120 °C in a fairly thin layer. (b) An initiation 
of downward induced flow from the bottom of central column toward the bottom of 
outer cylinder immediately creates a gradually expanding positive temperature 
gradient along the downward directions for time t>0.1 s. (c) A flow pattern of 
induced inflow from the surface of outer cylinder toward the center of axial 
coordinate promotes a parabolic shape temperature profile with downside tip point 
around the center of axial coordinate due to the induced downward cold streams 
from the bottom of central column for time t>0.1 s. (d) An induced upward 
profile of temperature distribution at the exterior region of container around 
the surface of outer cylinder is a result of induced upward flow pattern along 
the corresponding surface. (e) An induced downward elongation profile of 
temperature distribution immediately around the neighborhood of the surface of 
inner cylinder (central column) , is a result of induced downward flow pattern 
along the corresponding surface. 

(B) Convection Driven by Reverse Normal Direction (V>= 180°) Gravitational 

Acceleration 

Time animation of the development of induced flow fields associated with 
reverse normal direction (tJ>=180°) gravitational acceleration is illustrated in 
Figure 6 for convection induced flow vectors (in terms of directions and 
magnitude of induced flow velocities) in r-z plane. The maximum velocity for 
each subfigure is specified separately at the corresponding subfigure. This 
problem can be divided into two sections for flow fields induced by the reverse 
normal direction (tj>= 180°) gravitational acceleration. For z > 6 cm, the problem 
associated with horizontal heating from the outer walls of cylinder toward the 
central wall within a narrow cylindrical slit of space. For z < 6 cm, this is 
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a typical problem of vertical heating from the top (downward gravity direction). 
As the source material is located at the surface of z = 0, and the location of 
vapor deposition is designed at bottom horizontal surface of central column at 
z = 6 cm, the flow profile in z > 6 cm has very little effect on the flow pattern 
of vapor deposition. Thus, we can ignore the flow profiles induced within a 
narrow cylindrical slit of space. As to the flow in z < 6 cm, Figure 6 shows a 
stably stratified medium and natural convection is completely prohibited. 

For the purpose to illustrate the flow patterns of inflow and outflow 
induced by the reverse normal direction (V>=180°) gravitational acceleration. 
Figures 7 and 8 show time animation of initiation and development of outflow in 
the region below the central column (at z = 5.55 cm), and inflow in the region 
near the bottom of outer cylinder (at z = 0.45 cm), respectively. Figures 7 and 
8 also show flow vectors (in terms of directions and magnitudes of induced flow 
velocities) in r -6 (x-y) plane. The maximum velocity for each subfigure is 

specified separately at the corresponding subfigure. Comparison of Figures 6,7, 
and 8 show some interesting results as follows: (a) Practically, there is no 

induced flow field for natural convection along the vertical direction and only 
a stratified medium is shown in the r-z plane for z < 6 cm due to the phenomena 
of vertical heating from the top (downward gravity direction) . (b) There is a 
very small induced flow along horizontal direction in r-9 plane for z < 6 cm due 
to the heating from side walls, a slightly greater magnitude at top than that at 
bottom in this region for effect of low temperature central column at z > 6 cm. 

Figure 9 shows a time animation of the initiation and formation of 
temperature profiles driven by the reverse normal direction (ip= 180°) gravit- 
ational acceleration in r-z plane at 0=0° and 180° during the course of the 
entire flow region development. The curves show in Figure 9 are the isothermal 
lines with values indicated the temperature in °C at the corresponding locations. 
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Again, it shows a clear two regions. For region of z > 6 cm, horizontal heating 
from side wall drives the convective flow toward the central wall within a narrow 
cylindrical slit which shows a series of isothermal curves with temperature 

A 

decreasing from side wall toward the wall along central column. As indicated 
earlier, this region is not important in our study of vapor deposition from 
source materials located at the surface of z = 0 transferring to the bottom of 
central column at z = 6 cm. There is not much variation in the dynamics of 
temperature distribution for the region of z < 6 cm because the region is filled 
with stably stratified medium, and the natural convection is completely 
prohibited. 

(C) Convection Driven by 45° Oblique Direction Gravitational Acceleration 

In this case, 45° oblique gravitational acceleration is applied along the 
r-z plane with 0=0°. In other words, induced flow profile, shown in Figure 10 
with r-z plane at 0=0° and 180°, is an illustration of convective flow induced 
by a vector summation of gravitational acceleration with a magnitude of 
horizontal 0.707 g 0 (g Q sin 45°) rightward, and vertical 0.707 g 0 (g 0 cos 45°) 
downward components. On the otherhand, induced flow profile, shown in Figure 11 
with r-z plane at 0=90° and 270°, shows convective flow induced by a reduced 
gravity acceleration with a magnitude of 0.707 g Q downward component. 

Comparison of Figures 2 and 10 show the following results: (a) A downward 
flow toward the bottom of the outer cylinder with inflow stream is initiated at 
the bottom of central column for both normal and oblique gravitational acceler- 
ation at time t < 0.2 s. (b) A downward flow along the surface of cylindrical 
central column with inflow from the surface of outer cylinder toward the surface 
of cylindrical central column (at upper section of container within the slender 
gap between co-axial cylinders) is also initiated for both normal and oblique 
gravitational acceleration during the time period of 0,2 s < t < 0.4 s. (c) A 
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flow patterns of upward flow along the surface of left- side outer cylinder and 
downward flow along the surface of inner cylinder (central column) and right-side 
outer cylinder are gradually formed (at upper section of container within the 
slender gap between co-axial cylinders) with intensified flow velocities 
accompanied by a stream of inflow at left- side from the surface of outer cylinder 
to inner cylinder, and outflow at right-side from the surface of inner cylinder 
to outer cylinder driven by 45° oblique direction gravitational acceleration 
during the time period of 0.2 s < t < 0.9 s, This asymmetric flow profile is 
quite different from the symmetric flow profile resulted from direction 
normal gravity acceleration, (d) A tendency of inflow pattern at left-side, and 
outflow at right-side, for height in z > 1.2 cm; and outflow at left- side, and 
inflow at right-side, for height in z < 1.2 cm are gradually formed for time t 
> 0.9 s. (e) It shows a flow pattern of asymmetric recirculation flow with 
center at (r c , z c )=(0.6, 1.2) cm where r c and z c denotes the center of 
recirculation flow along r- and z- coordinates, respectively. Direction of 
recirculation flow is clock-wise for time t > 1.2 s. 

In contrast to Figure 10 for induced asymmetric flow profile in r-z plane 
at 0=O°-18O° due to ^=45° oblique direction gravitational acceleration, Figure 
11 shows an induced symmetric flow profile in r-z plane at 0=90°- 270° due to a 
reduced 0.707 g Q gravity acceleration in downward direction. Comparison of 
Figures 2 and 11 show the differences in induced flow fields driven by normal (1 
g Q ) and reduced (0.707 g Q ) gravitational accelerations. Major differences are 
as follows: (a) Upward motions along the surface of outer cylinder are shown for 
flow induced by the normal gravity acceleration, while the similar upward motions 
around the central axis of the cylinder are resulted for flow induced by the 
reduced gravity acceleration. (b) Downward motions around the central axis of 
cylinder are illustrated for flow induced by normal gravity acceleration, while 
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the similar downward motion along the surface of outer cylinder are resulted for 
flow induced by the reduced gravity. (c) There are two recirculation flows 
generated with centers at (r c , z c ) = (±1.9, 2.7) cm, and at (r c , z c ) = (±0.8, 3.2} 
cm for flows induced by normal and reduced gravity accelerations, respectively, 
(d) Directions of left-hand side and right-hand side recirculation flows are 
clockwise and counter clockwise, respectively, for flows induced by normal 
gravity, while that of the directions of recirculation flows for the similar 
sides are counter clockwise and clockwise, respectively, for flows induced by 
reduced gravity acceleration. These results show a very important result that 
the reduced gravity can help in promoting a direct transport of source materials 
at the bottom of outer cylinder to the location of vapor deposition at the bottom 
of central column without going through the complicated flow profiles of 
transport process . 

For the purpose to illustrate the flow patterns of inflow and outflow, 
Figures 12 and 13 show time animation of initiation and development of mostly 
outflow in the region right below the central column (at z — 5.55 cm), and mostly 
inflow in the region near the bottom of outer cylinder (at z — 0.45 cm), 
respectively. Figures 12 and 13 also show flow vectors (in terms of directions 
and magnitudes of induced flow velocities) in r-0 (x-y) plane. The maximum 

velocity for each subfigure is specified separately at the corresponding 
subfigure. It is interesting to emphasize that there is a flow region with 
mostly outflow for flow from the central axis of cylinder outward toward the 
direction of the surface of outer cylinder with location at the upper section of 
container with heights in z > 0.8 cm. It is also a flow region with mostly 
inflow for flow from the surface of outer cylinder inward toward the central axis 
of cylinder at the lower section of container with heights in z < 0.8 cm. 
Because horizontal component with a magnitude of 0.77 g Q gravitational 
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acceleration contributes to the initiation of very complicated flow profiles, 
shown in Figures 12 and 13, there are inflow component shown in height z > 0.8 
cm, and outflow component shown in height z < 0 . 8 cm in the range between 135°^ 

< $ < 225°. In other words, because of the presence of horizontal component 
gravitational acceleration due to oblique direction gravity acceleration, the 
vapor deposition at the bottom surface of the central column are contributed by 
the mixing of (a) inflow from external flow zones in the range between 135° < 9 

< 225° which were brought upward along the surface of outer cylinder from the 
evaporation of source materials stored at the bottom surface of outer cylinder, 
and (b) the direct transport from the source materials which were upward to the 
location of vapor deposition in the rest of the ranges. In other words, the flow 
analysis indicates that the vapor deposition at the bottom surface of central 
column was able to secure a complete result of direct straight upward diffusion 
from the evaporation of source materials stored at the bottom surface of outer 
cylinder exposed to oblique direction gravitational acceleration. In other 
words, obliqued direction gravitational acceleration can not secure a better 
quality vapor deposition because the presence of horizontal component of 
gravitational acceleration. 

Figures 14 and 15 show the time animation of the initiation and formation 
of temperature profiles in r-z plane at 0=O°-18O°, and at f?=90°-270°, respectively 
during the course of the entire flow region development under oblique direction 
gravitational acceleration. Figure 14 shows the temporal variation of isothermal 
lines resulted from the combination of rightward and downward gravity acceler- 
ations. It shows clearly the initiation and development of the progressing ox 
isothermal lines following the trend of a combination of rightward and downward 
movement. Comparison of Figure 10 for the velocity profile and Figure 14 for the 
temperature profile illustrate a very interesting similarity in these two groups 
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of curves. 


Figure 15 shows the similar temporal variation of isothermal lines from 
another angle point of view which was resulted from a reduced gravity 
acceleration with a magnitude of 0.707 g 0 . Initiation and development of cone 
shape temperature profile is resulted from a direct upward induced flow from 
bottom surface of outer cylinder transporting toward bottom surface of central 
column. These temperature profiles shown in Figure 15, and flow profile shown 
in Figure 11, suggest a very promising indication that the vapor deposition under 
a reduced gravity environment may be able to promote material processing of vapor 
deposition with a homogeneous crystalline films with fine grains. 

IV. Discussion and Conclusion 

In this study, mathematical formulations of the initiation and development 
of fluid profiles and temperature distributions driven by the heat transfer 
process of natural convection under normal and oblique directions of gravit- 
ational acceleration with extremely low pressure environment have been carried 
out. Purpose of present study is to investigate the characteristics of vapor 
deposition under low pressure environment to be used for promoting the process 
of nonlinear optical materials by growing thin crystalline films of diacetylene 
monomer (Nayer and Winter, 1990). Numerical computation of a full set of non- 
linear fluid dynamics show that vapor deposition process has to be accomplished 
through a fairly complicated flow patterns of recirculation under normal gravity 
environment both in normal and oblique directions which are the major driving 
forces to induce natural convection in low pressure environment. For the case 
of natural convection under normal gravitational acceleration with oblique 
direction along a plane perpendicular to the gravity acceleration, result shows 
a possibility to induce a simple direct updraft motion to transport source 
materials from the location of storage to the location of vapor deposition under 
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a reduced gravitational acceleration. Since the oblique direction gravitational 
acceleration is always accompanied with the component of gravitational 
acceleration in horizontal direction, the existence of horizontal component i$ 
responsible for the generation of complicated flow patterns with recirculation. 
It can be anticipated that there is no way to produce a homogeneous thin 
crystalline films with fine grains under such a complicated flow patterns of 
recirculation with a nonuniform temperature distributing in normal gravity 
environment. This study suggests that it is necessary to adopt an unique 
environment with microgravity to eliminate natural convection. In so doing, we 
can carry out an experiment which can simplify the flow pattern and reach a goal 
to produce a homogeneous thin films with fine grains in the future space 
experiments. Numerical study of fluid mechanics computation can be used to 
facilitate most optimistic way of experiment with best setup to achieve the goal 
of processing best quality nonlinear optical materials in the future space 
experiments . 
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Figure Captions 


Figure 1 


Figure 2 


Figure 3 


Figure 4 


Figure 5 


Figure 6 


Figure 7 


Figure 8 


Figure 9 


Figure 10 


Figure 11 


Figure 12 


Figure 13 


Geometrical illustration of co-axial cylinders to be used for 
mathematical modeling, (A) in r-z plane, and (B) in r-0 plane. » 
Flow profile initiation and development in r-z plane at 0=0 ° and 
180° under normal direction (V>=0 o ) gravitational acceleration. 

Flow profile initiation and development in r-0 plane at z = 5.55 cm 
under normal direction (V>=0°) gravitational acceleration. 

Flow profile initiation and development in r-0 plane at z = 0.45 cm 
under normal direction (^>= 0°) gravitational acceleration. 

Initiation and formation of temperature profile in r-z plane at 0=0° 
and 180° under normal direction (V>=0°) gravitational acceleration. 
Flow profile initiation and development in r-z plane at 0=0° and 180° 
under reverse normal direction (^>=180°) gravitational acceleration. 
Flow profile initiation and development in r-0 plane at z = 5.55 cm 
under reverse normal direction (^>=180°) gravitational acceleration. 
Flow profile initiation and development in r-0 plane at z = 0.45 cm 
under reverse normal direction (^lSO 0 ) gravitational acceleration. 
Initiation and formation of temperature profile in r-z plane at 0=0° 
and 180° under reverse normal direction (^lSO 0 ) gravitational accel- 
eration. 

Flow profile initiation and development in r-z plane at 0—0° and 180° 
under oblique direction (t/>= 45 0 ) gravitational acceleration. 

Flow profile initiation and development in r-z plane at 0=90° and 
270° under oblique direction (V>= 45°) gravitational acceleration. 
Flow profile initiation and development in r-0 plane at z = 5.55 cm 
under oblique direction 45°) gravitational acceleration. 

Flow profile initiation and development in r-0 plane at z = 0.45 cm 
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under oblique direction (V>= 45°) gravitational acceleration. 

Figure 14 Initiation and formation of temperature profile in r-z plane at 0=0° 
and 180° under oblique direction (V>= 45°) gravitational acceleration, 
Figure 15 Initiation and formation of temperature profile in r-z plane at 0=90° 
and 270° under oblique direction (V>=45°) gravitational acceleration 
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Flow Profile Initiation and Development in r - 0 Plane at z - 0.4 cm, \j/ - 



Initiation and Formatiom of Temperature Profile in r - z Plane at 0 = oo - 1800, y = 450 





Initiation and Formatiom of Temperature Profile in r - z Plane at 0 - 0° - 180° 









